Substrate temperature changes during MBE growth of GaMnAs 
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Remarkably big increase of the substrate temperature during the low-temperature MBE growth 
of GaMnAs layers is observed by means of band gap spectroscopy. It is explained and simulated in 
terms of changes in the absorption/emission characteristics of the growing layer. Options for the 
temperature variation damping are discussed. 
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I. INTRODUCTION 

The growth temperature is an important parameter de- 
termining the quality of the MBE grown ferromagnetic 
GaMnAs layers. Its optimum value depends on the Mn 
doping leveljii^i^ With increasing Mn density the growth 
temperature has to be reduced in order to prevent the 
Mn precipitation. At the same time, however, a certain 
minimum temperature must not be underrun in order to 
maintain the 2D growth. The growth temperature ad- 
justment is critical especially at high Mn doping levels 
where the temperature range of the 2D growth becomes 
narrow. Unfortunately, determining the substrate tem- 
perature at this temperature region is biased by a con- 
siderable uncertainty: the GaAs substrate is transparent 
for a low-temperature pyrometry, and a thermocouple 
can only weakly be linked to a rotating substrate holder. 
Moreover, it is difficult to eliminate a substantial delay 
in the temperature regulation. 

AppHcation of the diffuse reflectance spectroscopy^ (or 
band edge spectroscopy, BES) not only allows one to 
overcome the above problems, but also reveals remark- 
able temperature changes connected with the massive Mn 
doping. Similar observation has already been reported 
in connection with growth of smaller band gap layers on 
larger band gap substrates)^ and with heat flux transients 
due to opening/closing the effusion cell shutters^ The 
latter problem is certainly involved in our case, too. How- 
ever, as we show further, yet another mechanism takes 
effect, becoming essential when growing heavily doped 
samples mounted in In- free sample holders. 



II. EXPERIMENT 

The growth experiments were performed in Veeco 
Gen II MBE system. A semiinsulating SOO/um thick 
GaAs substrate with 2-inch diameter was mounted into 
a molybdenum In-free sample holder, fixed mechanically 
at its edge. No diffuser plate was used between the sub- 
strate and the heater. Distance between the heater (2- 
inch PBN plate with tungsten filament inside) and the 
substrate was about 10 mm. The band-edge spectrome- 
ter (kSA BandiT) was mounted on the central pyrometer 
port, normal to the substrate. In this arrangement the 
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FIG. 1: Temporal evolution of substrate temperatures mea- 
sured by the band-edge spectroscopy. The grown layers differ 
by the Mn contents of 2.5% (curve a), 3.5% (curve b), 5% and 
7% (curves c and d), respectively. Times f = (opening the 
Ga source), t = 60 s (opening the Mn source), and t — 1950 s 
(closing Ga and Mn sources) are indicated by vertical lines. 



thermal radiation of the heater had sufficient intensity 
near the band-gap wavelengths to serve as a radiation 
source for the transmission measurement even at tem- 
peratures as low as 200° C. 

The low-temperature part of the structures under in- 
vestigation consisted of 5 nm thick GaAs buffer and 
150 nm thick Gai_xMnxAs layer, with x varying from 
2.5% to 7%. At the growth rate used the temperature 
of the Ga source was 930°C, the Mn source tempera- 
ture varied from 770 to 830° C. During the growth the 
substrate heater was supplied by a constant power, with- 
out any feedback temperature control. The same initial 
substrate temperature of 182 ±2° C was reproduced when 
applying the same power and using always the same piece 
of sample holder. 

In Figure [1] temporal evolutions of the BES-measured 
substrate temperatures (Tbes) ^re shown for growths 
with various Mn densities. All the curves have several 
common features: (i) small (less than 5 degrees) increase 
of Tbes after opening the Ga source; (ii) pronounced (40- 
60 degrees) Tbes increase after opening the Mn source; 
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FIG. 2: Near band gap transmission spectra of an undoped 
GaAs substrate (curve a), and the same substrate with a 
50 nm thick GaMnAs layer with 7% Mn (curve b). Both 
curves were measured at room temperature. For contrast, 
curve (b) is repeated in the dashed line with a wavelength 
shift of AA = 19 /im which would correspond to increase of 
temperature from 200 to 250 °C (Ref.d). 



(iii) Tbes decrease after closing the sources, with a new 
steady-state value well above the initial substrate tem- 
perature. During all the growths the temperature mea- 
sured by usual floating thermocouple behind the sub- 
strate showed only a weak increase (less than 3 degrees) 
which relaxed back after the growths. 

When considering the dramatic difference between the 
Tbes response to opening the Ga and Mn sources, a ques- 
tion arises on the validity of the BES interpretation as 
a change in substrate temperature. An alternative ex- 
planation could consist in a change of the transmission 
spectrum shape due to presence of Mn impurity levels 
close to the fundamental absorption edge. If it were so, a 
relevant difference would have to be detectable between 
spectra of an undoped GaAs and a GaMnAs layer at the 
same temperature. Two such spectra are shown in Fig. [21 
representing a semiinsulating GaAs substrate, and the 
same substrate with a GaMnAs layer on top of it. In the 
band edge region, however, there appears to be no signif- 
icant difference between the spectra, comparable to the 
expected shift due to difference in temperature.—"^ Thus, 
the chemically induced spectrum change is unlikely in 
the relevant parameter range and the observed band edge 
shift can indeed be attributed to the change in tempera- 
ture. The more the question recurs on the mechanism of 
the dramatic temperature increase after opening the Mn 
source. 

We attribute this effect to increasing heat absorption 
in the heavily doped epitaxial layer. Two straightfor- 
ward absorption channels may be considered, related to 
the Mn doping: excitation of valence electrons to the 
Mn acceptor level, and free carrier absorption. The first 
channel can readily be neglected: with the ionization en- 
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FIG. 3: Infrared transmission spectra of undoped GaAs sub- 
strate and GaMnAs layers grown on the same substrate. Pa- 
rameters of the samples are as follows: (a) undoped substrate; 
(b) 50 nm, 7% Mn; (c) 150 nm, 7% Mn. AU curves were mea- 
sured at room temperature. 



ergy around 100 meV (Ref. d) the Mn atoms are almost 
fully occupied at the room temperature, making further 
transitions from the valence band impossible. On the 
other hand, the absorption efficiency via free holes is 
strong due to their extremely high density (in the or- 
der of lO^^cni^^). Since, moreover, their mobility is low 
and relaxation time short (in the order of 10 cm^/Vs, 
and 10~^^s, respectively), a significant absorption can be 
expected in a broad infrared range, starting already in 
the GaAs bandgap.^ This is illustrated in Fig. [3] where 
transmission spectra between 1 and 25 im\ are shown for 
various samples. 

Besides the free hole absorption the incoming heat is 
absorbed also by crystal lattice vibrations, independently 
of the doping. A strong two-phonon absorption is seen 
in the long-wavelength part of Fig. [31 This absorption 
and its one-phonon analog^" around 35 /im are mainly 
responsible for heating of an insulating GaAs substrateJ^ 



III. MODEL 

It is tempting to describe the heat exchange in the sub- 
strate simply as difference of (i) a heating term (propor- 
tional to heat flux from thermal sources, with an absorp- 
tion factor increasing in time), and (ii) a cooling term 
(proportional to deviation from a steady state tempera- 
ture). However, such a model turns out to be inadequate, 
yielding some results contradictory to experiments. For 
example, it predicts an unbounded increase of the sub- 
strate temperature, as long as the absorption increases, 
i.e. until the layer becomes completely opaque. Figures[Tl 
and [21 show a different behavior: the temperature stabi- 
lizes in spite of the nonsaturated absorption. The reason 
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of this failure is fundamental: according to Kirchhoff 's 
law the absorption efficiency of a body changes identi- 
cally to its emission efficiency at each wavelength,— i.e. 
no temperature change can occur unless the radiation 
spectrum changes. Disregarding the spectral distribution 
in a system with significant radiative coupling, as in the 
above model, inevitably leads to nonphysical solutions. 

A mathematical model we propose assumes fully radia- 
tive heat exchange. Its essential ingredients are proper 
description of the spectral power distributions of the ther- 
mal sources involved, and the doping induced change in 
the absorption and emission characteristics of the sub- 
strate. 

A radiating and absorbing body can be approximated 
by a black body with spectral radiance (monochromatic 
emissive power) i?(A) enclosed by a surface with dimen- 
sionless emittance e(A) and absorptance a(A). The radi- 
ance R{X) obeys the Planck law. 
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(1) 



The net heat transferred from a radiating body 1 to 
an absorbing/radiating body 2 can then be described as 

/•oo 

Qnet = / [^2 (A)ei (A)i?i (A) - £3 ( A)i?2 ( A)] rfA 

Jo 

(2) 

where dimensionless factor Fi_2 represents the geomet- 
ric configuration of the two bodies, 5*2 is the surface area 
of body 2. This general principle can be used to de- 
scribe — pair by pair — the heat exchange between the 
substrate and all other relevant thermal sources in the 
growth chamber. 

To simplify the description we adopt the following as- 
sumptions: 

(?) emittance of each body equals to its absorptance, 
e(A) = a(A) (diffuse form of Kirchhoff 's lawii); 
{ii) the heater and the cells behave like gray bodies, 
i.e. their emittances are constants; 

[iii) absorptance/emittance as (A) of the substrate is 
composed of the phonon and the free carrier contribu- 
tions. For simplicity, the phonons are assumed to absorb 
perfectly in their wavelength interval(s), (Api,Ap2)j the 
free carrier contribution a/ is assumed to rise in time, 
starting from zero in the undoped substrate and expo- 
nentially saturating as the total number of free holes in- 
creases, 



Ar>l < A < Ar) 



■pi ^ y\ ^ Ap2 /o^ 

af — 1 — exp(— fc/t), A elsewhere 



where the rate constant kf is proportional to dopant flux. 

The following bodies are relevant for the heat ex- 
change: the substrate heater (temperature T/;), the ef- 
fusion cells opened during the growth (effective temper- 
ature Tc), the background of the growth chamber (ef- 
fective temperature T;,), and the substrate (temperature 
Tg). The view factor F for each pair of bodies can be 



approximated by the solid angles A^,, A^, Af,, at which 
the substrate with area Sg "sees" the respective bod- 
ies; identity Ah = Ai, + Ac holds, as the heater covers 
one complete hemisphere of the substrate's scope (Ah), 
the other hemisphere being filled by the background (Ah) 
with small hot spots of the opened cells (Ac). The heat 
balance in the substrate is then 

Qnet ^ SsE A,[R{\; T,) - i?(A; T,)] a, dX, 

i 

i = h,b,c (4) 

Absorptance as (A) is a piecewise constant function, with 
its phonon part lying in a narrow interval in the mid in- 
frared range. There the Planck's law approximates to 
R oc T/X"^ (Rayleigh- Jeans law), leading to linear terms 
in T when integrating Eq. ([3]), whereas the complete 
Planck's law yields terms proportional to T^. Thence, 
equation ^ after integration results in 



Qnet — '^Ss ^ 

{af [T^ + a.T^ + (1 - ac)T^ - 2Tf] + 
(1 - a/)/3 [Th + acTc + (1 - ac)Tt - 2rs]} 



(5) 

where a is the Stefan-Boltzmann constant, = Ac/ Ah 
is the view angle ratio of the opened cells, and A},/ Ah = 
1 — flc. The constant factor (3 reflects the ratio of power 
irradiated in the phonon domain and power over all wave- 
lengths. 



{i/x% - i/xl,) 



(6) 



Taking Api = 16.5 /im and Ap2 — 18.5 /im leads to /3 = 
1 X 10^ K'' as an order of magnitude estimate,^'* 

In equilibrium, Qnet — 0, and Eq. ([5]) determines the 
steady state substrate temperature Tg for any given T^, 
Th, Tc, and a/; parameters (3 and ac are physical and 
technical constants, respectively. Out of equilibrium, the 
temporal evolution of is governed by 



CdTs/dt = Qnet, 



(7) 



where C is constant heat capacity of the GaAs substrate; 
varying electronic specific heat can safely be neglected 
at the relevant temperaturesi^^ Let us note that if the 
change in temperature is slow enough the system remains 
quasistationary and Tg{t) can still be find alone by solv- 
ing Qnet — with a/(t) increasing in time according to 
Eq. ©. 

Explicit solutions can be obtained in two special sta- 
tionary cases: 

a) Negligible free electron absorption, ov af — Q. Equa- 
tion ([5]) simplifies to 



Tg - {Th + acTc + (1 - ac)n)/2 



(8) 



In particular, Tg = {Th + T^) /2 in the steady state before 
the growth, when ac = 0. 
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FIG. 4: Calculated temperature evolutions of substrate 
temperatures of strongly and weakly doped layers (solid 
lines). Also shown are corresponding quasistationary solu- 
tions (dashed lines). 



b) Prevailing free electron absorption, or afTf ^ (1 — 
af)(5Ts. Then, 



r4=(T4 + a,T4 + (l-a,)r4) /2 



(9) 



In particular, = (T^ + T^)/2 in the steady state after 
the growth, explaining the persistent increase of Ts as a 
consequence of change of the absorption characteristic in 
presence of two different thermal sources. 

The above stationary solutions can be used to deter- 
mine the model parameters from the measured data. In 
accordance with Fig. [1] we can take Tg = 183°C before 
the growth, Tg = 243° C at the end of the growth of the 
most strongly doped layer, and Tg = 233° C long after 
closing the cells. Taking further = 0.005 as a rough 
estimate of the view angle ratio of the opened cells, we 
can solve Eqs. (0) and ®, getting Tt = 322°C for the 
substrate heater temperature, Tc — 954°C for the effec- 
tive temperature of the opened cells, and T;, = 44°C for 
the effective background temperature. This single set of 
plausible model parameters yields quasistationary solu- 
tions for Ts shown in Fig. [4] by dashed lines, if we chose 
A:y = 10 X lO^^s^i (2 X lO^^^s^^) as the doping rate 
constant of the strongly (weakly) doped layer. Trying, 
finally, aSsC-^ = 1 x 10"i°K-3s-i (corresponding to 2- 
inch GaAs substrate of 500 /xm thickness), Eqs. (O and 
^ lead to time evolutions Ts{t) plotted by solid lines in 
Fig.H 



IV. DISCUSSION 

On account of its simplifications the model cannot have 
the ambition of a perfect agreement. Its main virtue lies 
in the fact that it links up parameters of a straightfor- 
ward physical meaning with the observed temperature 



evolution, and that physically reasonable values of the 
model parameters lead to expectable results. 

The accuracy of the agreement could be arbitrarily in- 
creased, e.g., by including additional thermal sources in 
the growth chamber and by tuning their temperatures. 
This, however, would not necessarily make the model 
more realistic. There are numerous influences out of 
reach of a simple description. For example, a part of the 
heat exchange is due to a conduction between the sam- 
ple holder and other massive parts of the holder mech- 
anism, all of them with particular heat capacities; the 
free electron absorption is wavelength dependent in the 
IR range; all surfaces exposed to material fluxes change 
their absorptance/emissivity during the growth, leading 
to changes of their temperatures; neither the phonon do- 
main is absolutely absorbing, nor the near infrared ab- 
sorption is completely absent in the undoped substrate; 
reflection takes place, etc. 

The increase in the substrate temperature is hardly vis- 
ible for the floating thermocouple behind the substrate, 
and a temperature control based on it is thus inefRcient. 
Yet, it is desirable to keep the temperature change small, 
especially in highly doped GaMnAs layers, where the 
growth temperature is strongly constrained both from 
the bottom (disordered condensation) and the top (pre- 
cipitation). Several schemes can be applied to meet this 
goal: 

1. BES-locked substrate heater control; 

2. anticipatory throttling of the substrate heater; 

3. In-bonding of the substrate to a massive sample holder 
with inherently strong IR absorption; 

4. incorporation of a less doped (uncritical in temper- 
ature) but sufficiently absorbing GaMnAs sublayer be- 
tween the GaAs substrate and the final GaMnAs struc- 
ture. 

Temperature evolutions of the In-bonded substrate 
and the structure with the buried GaMnAs heating layer 
are shown in Fig. [S] The curves clearly differ in their 
time constants, corresponding to different heat capaci- 
ties. The temperature increase, however, is almost iden- 
tical in both systems, reduced to less than one half if 
compared with Fig. [TJ 

This illustrates an important difference between sys- 
tems possessing a significant IR absorptance (such as In- 
bonded substrates, or sandwiched structures with buried 
heating layer), and systems where a wide-range IR ab- 
sorptance is missing (such as In-free mounted semiin- 
sulating substrates). In the first case the temperature 
increase is solely due to activation of an additional ther- 
mal source; wide-range IR absorption is already prevail- 
ing and its further change (if any) does not affect the 
heat-exchange balance. In terms of the presented model, 
terms linear in T are negligible against terms T"*. 

In the second case the change of the spectrum of the 
substrate absorptance/emittance is essential; the change 
in temperature corresponds to a shift of the absorption 
maximum from longer to shorter wavelengths. In terms 
of the model, at first the main heat exchange takes places 
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FIG. 5; Substrate temperature evolution during growth of 
SMn = 7% layers on (a) the In-free mounted substrate, (b) the 
In-free mounted sandwich structure with a buried GaMnAs 
heating layer, (c) the In-bonded substrate. The buried heat- 
ing layer in (b) is 150 nm thick with a^Mn = 5%, separated by 
25 nm GaAs from the top GaMnAs layer. 



heavy Mn doping via the following scenario: 
{i) Before the growth, the undoped substrate absorbs the 
radiation from the thermal sources mainly in the mid- 
infrared range (phonon domain). 

[ii) The doping leads to strongly enhanced free-carrier 
absorption in the near-infrared range, shifting the de- 
cisive spectral absorptance towards the shorter wave- 
lengths. 

{Hi) More heat is absorbed from the high-temperature 
(short wavelength) sources. To reach a new heat- 
exchange equilibrium, the substrate radiance has to shift 
accordingly, i.e. its temperature increases. 

The proposed model, in contrast to a linear cooling 
model, explains all the main features of the observed ef- 
fect, while preserving the involved physical mechanisms 
disclosed. Using plausible values of the input parameters 
it leads to quantitatively correct results. It yields con- 
clusions relevant also for temperature stability during the 
MBE growth generally. Especially, it sheds a new light on 
behavior of In-bonded and In- free mounted substrates. 



through terms linear in T, later the T terms prevail. The 
corresponding temperature increase would occur even if 
no additional heat flux appears. On the other hand, in 
a system with negligible absorption towards the shorter 
wavelengths, opening a high temperature cell with small 
orifice has little or no effect on the absorbed heat. Thus, 
an interesting paradox can occur: temperature variations 
in an In-bonded substrate can turn bigger than those in 
an In-free mounted substrate with weak free-carrier ab- 
sorption. 

V. SUMMARY 

The main part of the dramatic increase in the substrate 
temperature could unambiguously be connected with the 
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